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An  Expansion  of  the  Gegenbauer  Polynomial  C^(xy) 
Introduction 

Sonine’s  second  finite  integral  [1,  p.  376]  may  be  written 

n/2 

Jj/x  sin  0)  ix(y  cos  0)  sin'fH  0  cosA+l  0  d0 

Wi  (v^+y2 ) 

(VWy2  y+A+1 


X 


for  ail  complex  x  and  y,  and  is  valid  when  both  Re(pi)>-1  and  Re(A)>-l.  At  least  two 
proofs  of  this  result  are  known.  One  involves  expanding  the  integral  in  powers  of  x 
and  y;  the  other  involves  integration  over  subsets  of  the  surface  of  the  unit  sphere  in 
R3.  Both  are  given  in  [1]. 

For  the  case  of  real  m  and  A,  a  third  proof  is  given  here  that  depends  in  an 
essential  way  on  the  identity  (7).  In  this  connection,  the  particular  form  of  the 
coefficients  ak  n(y)  is  important;  that  is,  the  easily  derived  identity  (10)  does  not 
seem  to  be  at  all  useful,  but  the  identity  (8)  is  exactly  what  is  needed.  It  facilitates  the 
investigation  of  the  limiting  form  (27)  of  ak  n(y)  as  n  tends  to  infinity.  The  identity 
(8)  is  apparently  new;  however,  the  special  case  of  y  =  1  was  known  to  Gegenbauer. 

Equation  (8)  is  interesting  in  another  regard  as  well.  A  simple  inspection 
suffices  to  prove  that  akn(y)>0  for  all  n  and  k  whenever  y>l  and  n^A>0.  The 
coefficients  remain  positive  in  the  two  limiting  cases  ^  >  0,  A  =  0  and  n  =  A  =  0,  as  can 
be  seen  from  (18)  through  (21).  In  fact,  it  was  only  this  positivity  result  that  the 
author  originally  sought. 

The  result  (3)  of  Mehler-Heine  type  is  apparently  new.  It  is  needed  to  prove  (I) 
by  our  methods.  It  has  additional  interest  in  that  it  duplicates  Szego’s  result  (2) 
simply  by  setting  y  =  0.  Since  Szego’s  proof  of  (2)  may  very  nearly  be  lifted  verbatim 
to  prove  (3),  it  is  perhaps  surprising  that  he  does  not  mention  (3)  in  [2].  The  special 
case  (4)  involving  Chebyshev  polynomials  is  particularly  striking  and  seems  to  be 
new  also. 
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Derivations  and  Results 


Let  a  and  j3  be  arbitrary  real  numbers.  For  any  complex  number  x,  the  Mehler- 
Heine  theorem  states  that 


lim  n“°  pi  (cos|-)  =  (x/2)'9  Ja(x),  (2) 

n— oo 

where  Ja(x)  is  the  Bessel  function  of  the  first  kind  of  order  a  [1,  §3.1(8);  2  (1.71.1)]. 
A  straightforward  proof  of  (2)  can  be  found  in  Szego  [2,  Theorem  8.1.1].  Szego’s 
proof  can  be  readily  modified  to  show  that 


lim  n~° 

n 

n-*°° 


cos  n  )  =  (>/2  \Zx2“-yT  )-«  J  (Vx2-y2  ) 

/ 


for  all  complex  x  and  y.  Like  the  Mehler-Heine  result,  this  formula  holds  uniformly 
for  x  and  y  in  every  bounded  region  of  the  complex  plane.  The  special  case 
a  =1 3  =  -1/2  gives  the  interesting  result 

lim  T  (  cosTr  \  =  cos  Vx2-v2  ,  (4) 


where  Tn(x)  is  the  Chebyshev  polynomial  of  the  first  kind  [2,  (4.1.7)].  This  follows 
from  (3)  by  using  Stirling’s  formula  and  the  well  known  result  [2,  (1 .71.2)] 


J-1/2(Z)  COSZ’  (5) 

We  will  need  another  special  case  of  the  general  result;  specifically,  for  m>-1  . 

lim  q  (  cosJ*\  =  v£72  , 

n~~  2 iu  "  \  cos-*;-  I  2^r(H+l)(Vx2-y2  TVl  (6) 

where  C£(x)  are  the  ultraspherical,  or  Gegenbauer,  polynomials  [1,  (4.7.1)].  (Szego 
uses  the  notation  pW(x)  instead  of  CJJ(x).) 

We  derive  Sonine’s  second  finite  integral  by  finding  an  alternate  form  for  the 
left-hand  side  of  (6).  This  requires  the  following  result.  For  n>k>0,  the  coefficients 
ak  n(y)  in  the  expansion 


C>y)  =2  aM(y)  ci_2k(*)>  n  =  0, 1 , 2,  .  .  . 


are  given  explicitly  by 
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ak.n(y)  =  (n-2k+A)  (|4)n_k 


£  Qi-A  +  m)k  m(y2-l)m  y"-2m 
m=°  m!(k-m)!(A)n_k_m  +  1 


(8) 


where  we  take  0°  =  1  and  (0)0  =  1  whenever  they  occur.  Setting  y  =  1  in  (8)  gives 

ak  n(l)  =  (n~2k  -H  A)  (M)n-k  t  (9) 

kUA)n.k  +  1 

which  is  due  to  Gegenbauer  [2,  (4.10.27)].  Furthermore,  for  real  y>  1  and  ji>A>0, 
the  coefficients  ak  n(y)  are  all  positive  as  can  be  seen  by  inspection  in  (8). 

The  formula  (8)  is  derived  as  follows.  Let  n>k>0.  In  the  expression  [2, 
(4.7.31)] 


[n/2] 

Cftx)  =  X  (-1)" 


m  =0 


ml  (n-2m)l 


(2x) 


n-2m 


we  replace  x  with  xy,  substitute 


£  n-2m  J 

(2x)n-2m  ^  (n-2m  +  A-2s)  . 

- =  4.  — — -  '“n-2m-2s'x'’ 

(n-2m)l  s=o  s!  (A)n-2m-s+ ! 


and  collect  terms  to  get 


n-2m 


H  „(y)  =  I  (.lrln-2k-H)0,)„.my 


=  yn'2k  (n-2k  +  A)  Q.  (2y2-l), 


where  Q.  is  a  polynomial  defined  for  general  complex  argument  u  by 


QA  u)  =  I 


H)m  C m), 


m=0  m!  (k-m)!  (A) 


^ 

'n-m-k  + 1  \  ^  / 


GO) 

(11) 


(12) 


For  arbitrary  a  and  0,  the  Jacobi  polynomial  of  degree  k  >  0  can  be  written 


P^(u)  =  2  (~l)m  (k4~ g  +  0  +  l)k,m (k-m  +  0+  l)n 


m  =0 


m!  (k-m)! 


/u+l\k-m 
{-)  ' 


(13) 


«p«p«w 


which  follows  from  [2,  (4.21.2)]  using  the  identity  [2, (4. 1.3)].  Setting  a  =  p-A-l  and 
p  =  A  +  n-2k  in  (13)  shows  that 


Qk(u)  = 


n-k  +  ] 


pG*-i-l,A+n-2k^u) 


(14) 


Expanding  the  Jacobi  polynomial  in  (14)  using  [2,  (4.3.2)] 

pM)(u)  =  i  d+*\(i+i»k  f^±\m  /u+iy-mt 

k  m!  (k-m)!  (1  +  a)m  (1  +P)k_m  \  2  /  \  2  ) 


(15) 


and  substituting  u  =  2y2-l  gives 


Qk(2y2-D  =  0i)n.k 


2  (M-A  +  m)k_m(y2-iry^-2m 

m=°  m!(k-m)!(A)n.k.m  +  i 


(16) 


Thus  (16)  and  (11)  establish  (8). 

Two  limiting  cases  of  (7)  are  easily  derived  from  [2,  (4.7.8)] 

Jim  Y\  CnA(x)  =  Tn(x)>  ”>1,  U7> 


and  are  worth  recording.  Thus,  for  p  >  0, 


fn/2]  ^ 

Cftxy)  =  J  b 

k  =  0 


k,„(y)Tn-2kW. 


n  =  0,  1, 2, 


where 


and 


•VnW  “  2M„-k  Z 


k 

E 

m  =0 


Oi  +  m)k-m(y2-l)myn-201 

m!  (k-m)!  (n-k-m)! 


[n/2] 


T (xy)  -  I  ck,„(y)  Tn_2k(x),  n=  1,  2,  3, 


k  =  0 


where 


ck,n(y)  =  n(n-k-l)! 


k 

I 

m  =  0 


(m)k-m(y2~1)m  yn'2m 

m!  (k-m)!  (n-k-m)! 


(18) 


(19) 


(20) 


(21) 
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The  notation  S'  means  that  1/2  the  last  term  in  the  sum  is  taken  if  n  is  even,  and  all 
of  it  is  taken  if  n  is  odd.  Note  that  inspection  shows  that  y  >  1  implies  bk  n(y)  and 
ck  n(y)  are  positive. 

Sonine’s  second  finite  integral  is  now  derived  from  (6).  Fix  x  and  y.  Let 
N  =  [n/2] .  From  (7) 


n1'2^  r* 


cos7T 


1  ^  An1-2*4  (1 +N) 

1+N  k=o  /i(n-2k)‘-2X 


I  \  (n-2k)'~2x  c^co^l 

iir  /  2A 


=  f  fn(l-t)gn(l-C)d(, 


where  we  have  defined  for  0  <  £  «  1 , 

£  >--_*■(■+ n> 

‘■°  „(n-2k)'-“ 


g„o-«)  =  i  ("~-2t)l-;A  c;.;„(cosi)xEk(i-{), 


and  XEk  is  the  characteristic  (indicator)  function  of  the  interval 


_L  ,  ,  k  =  0, 

N  + 1  N  +  1  / 

'JL  ,  lull  ,  k  =  N 
N+l  N  +  lJ 


»  k  =  0,  1 . N-l, 


It  can  be  verified  that  XE(,(2k/n)  =  1  for  k  =  0,  1,  .  .  .  ,  N. 

Assume  for  the  moment  that  both  ] f n(C) I  and  |gn(£)|  are  bounded  above  by 
integrable  functions  of  £.  To  do  this,  it  will  be  seen  that  we  must  restrict  attention  to 
A > -1/2,  m  > -1/2,  p>k,  so  that  the  integral  [2,  (1.7.4)] 

f  Vi-^-'dt  =  rfr+'/orpi-*)  (23) 

0  2f0i  +  '/:) 
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will  be  finite.  If  f  =  lim  f  and  g  =  lim  gn,  the  bounded  convergence  theorem  [3,  p. 
1 10]  implies 

iim  Ji!^_  CM  /  cosTT 
—  2,  "  ycos-^- 

Let  £in  (0,1)  be  rational.  Then  l-£  =  2k/n  for  sufficiently  large  k  and  n,  so  that 
g(l-0  =  Iim  gn(l-£) 


lim 

n-*oo 

l-£  =  2k/n 


(n— 2k)1  2X  „jt  ,  x  . 

Cn-2k(COSlT) 


2A 


=  V^t/2 


2*f(A  +  1) (£x) 


\-Vi 


(25) 


with  the  last  step  following  immediately  from  (6).  Thus,  (25)  holds  for  all  £  in  [0,1] 
by  continuity.  Similarly,  from  (8)  and  for  all  £  rational  in  (0,1), 


f(l-£>  =  Hm  fn(l-£) 

n—oo 


=  lim 

n— «> 
l-{  =  2k/n 


An'-^l+N) 

H(n-2k)N2A 


=  lim  2! 

m  =0 

n-»oo 

1  -{  =  2k/n 


(£ ■ +  I )  o* ■ +  1  w,0'2*1'  l-,: ) (1  +  -I?)  0^~A  +  m)k_m  sin2m £ 


£|-2ACOsn  m!  (k-m)!  (A  +  l)n  k 


(26) 


Interchange  the  limit  and  the  summation,  and  evaluate  the  limit  of  the  mlh  term 
(convert  Pochhammer  symbols  to  Gamma  functions,  apply  Stirling’s  formula,  and 
use  k(n-k)  =  ( 1  -£2)  n2/4)  to  obtain 


f(l-0  = 


I 


m  =  0 


£2*(l_£2)M-*-> 


22m-^a-i 


T(A+1)  (’/2y  \ZT-£^)2ni 

r'0-<  -t- 1 )  mirOi-A  +  m) 


£2a(i_£2)m-*-i  r(A+i) 

22m-2A-1  TOi+1) 


U-,(yVT7) 

C/iyx/r^7)^-1 


(27) 


where  Iv(z)  denotes  the  modified  Bessel  function  of  the  first  kind  of  order  v  (see  [1, 
§3.7(2)].  We  must  require  ji>A  in  (27)  to  have  convergence.  Continuity  again 
assures  that  (27)  holds  for  all  £  in  (0,1).  Now,  interchanging  the  limit  and  the  sum 
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was  valid  because  an  upper  bound  for  the  total  sum  can  be  found.  Since  the  absolute 
value  of  the  m,h  term  in  (26)  is  bounded  by 


B  r(A+1) 

t(m+d 


c^iyi^r" 

m!  Rpi-A  +  m) 


where 


B  = 


M-2/l 


n-2<,i-A-l) 

n2m|cosn 


f(k  +  pi-A) 
Rk-m  +  1) 


Rn-k  +  pi) 
Rn-k  +  A  +  1-m) 


n—°°. 


the  total  sum  in  (26)  is  bounded  by 

p(£)  =  L£2A(1-£2)M-A_I  +  ^  X  C^lyl  )  <  op 

Rpi  +  1)  m  =  0  m'.rOi-A+m) 


(28) 


for  some  constant  L  independent  of  £.  The  series  in  (28)  is  a  continuous  function  of 
4  on  [0,11  >f  pi>A.  Hence,  from  (23),  F(£)  is  an  integrate  function  that  bounds  |f(£)| 
for  all  n. 


From  (24),  (25),  and  (27)  we  have 

lim  C"  /cosM 

n^°°  ycosn“  ) 


=  _ _ _  j ,1/2(tx)i  A_, (yv/i^F ) d^ 

2^  r(pt  +  1 )  xA~  ‘/2  yM-A-l  0 

^Tl  J  ,/2(VxM^) 

2M  T(pi+  1)  (Vx2-yJ  f~ 1/2 

with  the  last  equation  from  (6).  Substituting  £  =  sin0  and  y  =  iy'  in  the  last  two 
formulas,  and  setting 


pi'  =  A-'/z  >  -1  and  A'  =  pi-A-1  >  -1  (29) 

yields  Sonine’s  second  finite  integral  (1).  The  only  thing  left  to  prove  is  that  |gn(£)l  is 
bounded  by  an  integrable  function  on  [0,11.  Szego’s  argument  [2,  p.  192]  in  the 
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proof  of  (2)  can  be  modified  easily  to  show  jgn($)|  is  bounded  by  a  constant. 


Conclusions 

The  special  case  p  =  X  in  (27)  may  merit  further  study.  In  the  more  restrictive 
case  p  =  A  =  0,  it  is  known  that  [4,  (87 1 .2)1 

;Jl/2 

1  ,(y  cos  0)  cos(x  sin  0)  d0. 
o 

However,  we  do  not  pursue  this  here. 

Another  question  that  we  do  not  investigate  here  is  the  expansion 
PW)(xy )  =  £  Ak  (y) 

n  k  =0 

It  would  seem  difficult  to  obtain  a  form  for  Ak  n(y)  from  which  it  is  directly  evident 
which  conditions  imply  Ak  n(y)  >  0. 

The  proof  of  (1)  presented  here  was  intentionally  restricted  to  real  p  and  A. 
However,  it  is  not  hard  to  see  from  (23)  and  (29)  that  the  proof  can  be  carried  out 
for  complex  p  and  A  provided  appropriate  remarks  are  made  in  appropriate  places 
about  the  complex  case.  If  such  remarks  are  made,  our  derivation  proves  (1)  for 
Re0i)>-1  and  Re(A)>-l.  Divergence  of  (23)  is  seen  to  be  the  cause  of  the  restric¬ 
tions  on  p  and  A. 
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NAVSEADET  NORFOLK  1 

NASC  1 

NAVAIRDEVCEN,  Warminster  2 

NOSC  2 

NAVWPNSCEN,  China  Lake  1 

NCSC  1 

CIVENGRLAB  1 

NAVSURFWPNCEN  1 

CHESNAVFACENGCOM  1 

NAVPGSCOL  1 

NAVWARCOL  1 

APL/UW,  SEATTLE  1 

ARL/PENN  STATE,  STATE  COLLEGE  1 

ARL/UNIV  OF  TEXAS  1 

DT1C  12 

DARPA  3 

NOAA/ERL  1 

NATIONAL  RESEARCH  COUNCIL  1 

WOODS  HOLE  OCEANOGRAPHIC  INSTITUTION  1 

ENGINEERING  SOCIETIES  LIBRARY, 

UNITED  ENGINEERING  CENTER  1 

MARINE  PHYSICAL  LAB,  SCRIPPS  1 
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